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Abstract 

We consider the gradient flow of a one-fiomogeneous functional, wliose dual involves the 
derivative of a constrained scalar function. We show in this case that the gradient flow is 
related to a weak, generalized formulation of a Hele-Shaw flow. The equivalence follows from 
a variational representation, which is a variant of well-known variational representations for 
the Hele-Shaw problem. As a consequence we get existence and uniqueness of a weak solution 
to the Hele-Shaw flow. We also obtain an explicit representation for the Total Variation flow 
in one dimension, and easily deduce basic qualitative properties, concerning in particular the 
"staircasing effect". 

1 Introduction 

This paper deals with the L^-gradient flow of the functional 

Jk{uj).= I \Mdx G {0,...,A^- 1} 

J A 

defined on differential forms uj e L'^{A,VI^{MJ^)), where A C is an open set. We will focus 
on the particular case fc — — 1: in that case, the dual variable is a scalar and this yields very 
particular properties of the functional and the associated flow. 

Notice that, when fc = 0, the functional Jq reduces to the usual total variation. When k = N— 1 
we can identify by duality w G r2^"i(M^)) with a vector field u G R^), so that Jn-i 

is equivalent to the functional 

V{u) / \diYu\dx (1) 

J A 

that is, the total mass of div u as a measure. 

The gradient flow of T) has interesting properties: we show in particular that it is equivalent to 
a constrained variational problem for a function w such that Aw = div u. Moreover, under some 
regularity assumption on the initial datum uq, such a variational problem allows to define a weak 
formulation of the Hele-Shaw flow [9l [11] (see also [T2| for a viscosity formulation). Therefore, it 
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turns out that the flow of ((T|) provides a (unique) global weak solution to the Hele-Shaw flow, for 
a suitable initial datum uq. But our formulation allows us to consider quite general initial data 
Mo, for which for instance divug may change sign, or be a measure. 

The plan of the paper is the following: in Section [2] we introduce the general functional we are 
interested in, we write the Euler-Lagrange equation for its Moreau-Yosida approximation and, in 
Section [231 we express it in a dual form that will be the base of our analysis. 

In Section [3] we focus on the case k — I which is analyzed in this paper. We show many 
interesting properties of the flow: comparison, equivalence with a weak Hele-Shaw flow if the 
initial datum is smooth enough, and qualitative behavior when the initial datum is not smooth. 
In Section im we observe that, in dimension 2, the case fc = — 1 also covers the flow of the L^- 
norm of the rotation of a vector field, which appears as a particular limit of the Ginzburg-Landau 
model (see [ini [T5] and references therein) . 

Another interesting consequence of our analysis is that it yields simple but original qualitative 
results on the solutions of the Total Variation flow in dimension one (see also [3l [5]). We show 
in Section 14.21 that the denoising of a noisy signal with this approach will, in general, almost 
surely produce a solution which is "flat" on a dense set. This undesirable artefact is the well- 
known "staircasing" effect of the Total Variation regularization and is the main drawback of this 
approach for signal or image reconstruction. 

2 Gradient flow 

Given an initial datum uiq e L'^{A, r2'^(R^)), the general theory of [6| guarantees the existence of 
a global weak solution ui e i^([0, -l-oo), L'^{A, il*'(M^))) of the gradient flow equation of Jk- 



where dJk denotes the subgradient of the convex functional Jk- Given e > and / G L'^{A, 57'^(M^)), 
we consider the minimum problem 



t e [0, +oo) , 



(2) 




(3) 



Notice that 




The Euler-Lagrange equation corresponding to ([3]) is 



e{f-io)edJk{io), 



that is there exists a (fc -|- l)-form v with \v\ = 1 such that v 



da;/|da;| if duj ^ 0, and 



£(/- w) = d*v in A 



and 



[*v)t = on dA. 



(4) 



2.1 Dual formulation 



Equation (|4]) is equivalent to 



iuedrkieif^u;)), 
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where 

TU \ f ^ T ( ^ I ^ if hll*<l 

Jk{m-= sup / -q ■ w dx ~ Ju{w) = < 

w.a^vl^ J a y +00 otherwise 

and 

||r/||* = sup |y rj-wdx: Jk{w) < I 
Note that 

Jk (w) + Jfe (??) > / w ■ rjdx 
J A 

for all w, rj. The equality holds iff J^r/ ■ w dx = Jfe(w), and in such case we have ||?7||* < 1. 
Letting m be a minimizer of ([3]) and t] = (f — u)/e we then get 

u ■ — dx = Jkiu)^ 

A £ 

which implies 

Jl(^l—^^^0 that is e>||/-M||*. 

In particular, we showed the following (see also [15] for the same result in the case of the Total 
Variation). 

Proposition 2.1. The function u — is a minimizer of ^ if and only if 

£>£c:= 11/11*. (5) 

Note that ||ry||» < 00 implies that 

rjw — 



for all w such that dw = 0. By Hodge decomposition, this implies that ij — d*g for some 2-form 
g, with gN = on dA. It follows that 

||?7||* = sup / d*g-wdx= sup / g-dwdx+ / wA*gN= sup / g-dwdx. (6) 

/^|dtu|<lJyl lA\dw\<lJA J OA /^|du)|<l"'A 

We then get 

h\\* = inf llgllL-(A)- 

sn\oa=o 

Indeed, it is immediate to show the < inequality. On the other hand, by Hahn-Banach Theorem, 
there exists a form g' , with d*g' = d*g = 77 such that 



sup / g-dwdx— sup / g' ■ ip dx = \\g\\L^(^Ay 



f^\dw\<lJA /^I^I<1 

Fix now (j)Q such that d*(f)Q = We can write g — ipo + d*ip, so that (|6]) becomes 

WvW*^ min 1100 +d*V'||Loo(A)- (7) 
The Euler-Lagrange equation of ([7]) is similar to the infinity laplacian equation 

doo('/>o + d*V) - 0. 

By duality problem ([7]) becomes 

min ||VV' + 0o||l~(^), (8) 
and the corresponding Euler-Lagrange equation is 

((v^v^ + V0o)(vv + (vv + <^o)) = 0. (9) 



3 



3 The case k = N -1 



In this case, we recall that we are considering the gradient flow of the functional ([T]), which is 
defined, for any u e Ll^^{A;R^), as follows 



V{u) = sup -uVw da; : v e (A) ,\v{x)\ < 1 \/x e 



(10) 



This is finite if and only if the distribution divu is a bounded Radon measure in A. We now see 
it as a (convex, l.s.c, with values in [0,+oo]) functional over the Hilbert space L^(A;]R^): it is 
then clear from (fTO|) that it is the support function of 

K = {-\7v : z;ei7i(A;[-l,l])} 

and in particular p G d'D[u), the subgradient of T) at u, if and only \i p £ K and Jj^p ■ udx — 
V{u)=J^\divu\: 

dV{u) = |-Vt> : V e Hl{A:[-l,l]) , J -Vvudx = J |divu|| . 
We can define, for u e domP, the Radon-Nikodym density 

QIV U , , ,. JB(x,p) 

Odivu{x) = — -Ax) = hm 



|divw| ''^o /i3(x,p) Idiv^^r 

which exists |divw|-a.e. (we consider that it is defined only when the limit exists and is in {—1, 1}), 
and is such that divu = 0div«|divu|. We can also introduce the Borel sets 

= {x(^A: 0div«(2:) -±1} . 

Then, we have: 
Lemma 3.1. 

dV{u) = {-Vv : veHl{A- [-1,1]), v = ±l |divw|- a.e. on f ±} . 

Proof. Consider v G Hq{A; [—1, 1]). Then we know [T] that it is the limit of smooth functions 
Vn G C^{A; [—1, 1]) with compact support which converge to v quasi-everywhere (that is, up to a 
set of i/^-capacity zero). 

We recall that when u £ L'^{A;M.^), the measure divu G H~^{A) must vanish on sets of H^- 
capacity [2 §7.6.1]: it follows that jdivuj-a.e. in A. Hence, by Lebesgue's convergence 

theorem, 

\7v{x)u{x)dx = lim / Vn{x)6divu{x)\divu\{x) ~ / t;(a;)div u(a;) . 

It easily follows that if u = ±1 |divu|- a.e. on — Vu G dT>{u) and conversely, that if u G d'D{u) 
then V — ±1 |divu|-a.e. on £^ . □ 



We now define, provided u G douydV (i.e., d'D{u) ^ 0), 

d°'D{u) = argminjy \p\'^ dx : p £ d'D{u)'^ : 
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it corresponds to the eleraent p = — Vw G d'D{u) of minimal i^-norm. Using Lemma l3.ll equiv- 
alently, v is the function which minimizes jVwpiix among ah v £ Hq{A) with v > X£+ and 
V < — Xf-i |divu|-a.e.: in particular, we deduce that it is harmonic in A \ £u U £u ■ 
Let us now return to the flow ([2]) . In this setting, it becomes 

u{0) = uo 



where v satisfies bl < 1 and 



T){u) + / u-Vv = 0. 
J A 



It is well know, in fact, that the solution of (fTTI) is unique and that —Wv(t) ~ d'^'D{u{t)) is the 
right-derivative of u{t) at any i > [Q. Given the solution {u{t),v{t)) of (fTTj) . we let 

wit) := / v{s) ds, 
Jo 

which takes its values in [~t,t]. We have 

u{t) = Uo + Vw(t). 

Theorem 1. Assume uq € L'^{A;R^). The function w{t) solves the following obstacle problem 



^ J \uo + Vwp dx : w & Hq{A) ,\w\ <t a.e.j 



(12) 



Observe that in case we additionally have divwo > a > 0, this obstacle problem is well-known 
for being an equivalent formulation of the Hele-Shaw problem, see [51 lllj. 

Proof. Given uq G L'^{A;M.'^), we can recursively define G i^(^;R^) as the unique solution 
of the minimum problem 

min Vs{u,Un), 



Ve{u,v) ^V{u) + I -^\u-v\'^dx. 



where 

Then, there exists G 9I?(u„+i) such that 

Un+l -Un- eVVn+l = 0. (13) 

It follows that Vn+i G Hq{A) minimizes the functional 

2 



\un + eVvj dx 

J A 

under the constraint \v\ < 1. Let now 

n 

Wn := ey^^Vj 

i=l 

The from (fT3|) we get 



Un=Un+VWn, (14) 
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and Wn minimizes the functional 

/ luo + Vwpdx (15) 
Ja 

under the constraint \w — w„_i| < s. Notice that — Wn-i \ < £ for all n imphes 

\wn\ < ne. (16) 

We now show that Wn minimizes also under the weaker constraint ((T5)) . Indeed, letting Wn 
be the minimizer of (jl5[) under the constraint (jl6p . we have 



Wn~ £ < lin+l < Wn + £, 

which follows by noticing that min{w„, + e} and maxjwn, t«„+i — e} minimize (|15p. hence 
they are both equal to w„. It then follows w„ = w)„ for all n. 

Passing to the limit in n we get the corresponding result in the continuum case. □ 

Remark 3.2. The previous proof also shows that for any initial uq G L'^{A;R^), u{t) — uq + 
Vw(t) is the unique minimizer of 

|divu| "I" ^1'" ~ ""op dx. 

We recall that obviously, such property does not hold for general semigroups generated by the 
gradient flow of a convex function. It is shown in |2J to be the case for the Total Variation flow, 
in any dimension, when the initial function is the characteristic of a convex set. 



3.1 Some properties of the solution 

A first observation is that t w(t) is continuous (in Hq{A), strong), as follows both from the study 
of the varying problems (IT2|) and from the fact that the flow u{t) — uq + Ww{t) is both continuous 
at zero and i^(yl)-Lipschitz continuous away from t = (and up to f = if mq G dome*!?). 

In fact, one can check that w is also i°°-Lipschitz continuous in time: indeed, it follows from 
the comparison principle that for any s < t, 

w{s)-t + s < w{t) < w{s)+t-s (17) 

a.e. in A, hence |jz«(i) — u'(s)||^=o(^') < |i — s|. The comparison (jl7p is obtained by adding the 
energy in ()12|) of mm{w{t), w{s) + t — s} (which is admissible at time t and hence should have an 
energy larger than the energy of w{t)) to the energy of max{w{t) — t + s, w{s)} (which is admissible 
at time s), and checking that this sum is equal to the energy at time t plus the energy at time s. 
This is quite standard, see [71 [T2]. 

In particular, we can define for any t the sets 

E+{t) = {w{t)^t} and E-{t) ^ {w{t) ^ -t} , (18) 

where 'w{t) is the precise representative of w{t) G H^{A), defined quasi-everywhere by 

w{t,x) = lini — / w{t,y)dy (19) 
P^o ujnp'^ Jb(x,p) 



6 



(wat is the volume of the unit ball). It follows from ((T7)) and ([T5)) that if w{t, x) = t, then for any 
s < t, X is also a point where w^s^x) is well-defined, and its value is s; similarly if w{t,x) ~ —t 
then ■u;(s, x) = —s. Hence: the functions 1 1-^ E^{t), 1 E~{t) are nonincreasing. 
Also, if s < one has from (|T7)) 

— / w{s,y)dy ~t + s < — ^/ w{t,y)dy < — ^/ w{s,y)dy - s + t 

(^NP" Jb{x,p) f^NP" Jb{x,p) ^NP" Jb{x,p) 

SO that if X e E^{s), 

..If . 1 /" 
2s — t < liminf -rr / w{t,y)dy < limsup jrr / w[t,y)dy < t 

P^O CJArp^^ yS(2;,p) ^^WP^^ J B{x,p) 

and sending s to i, we find that if x e {~\s<t-^^ i^)^ w{t,x) — t and x E E^{t): hence these sets 
(as well as E^{-)) arc left-continuous. 
We define 

E+{t) = \Je+{s) CE+{t) and E' {t) - \J E' {s) C E' {t) , (20) 

s>i s>t 

as well as £■(*) = i;+(t) U ^^"(i), ^'^(i) = E+{t) U £'7(t). Then, there holds the following lemma: 
Lemma 3.3. If s <t, then 

E~{t) C E- (s) and E+ (t) CE+{s), 

E-{t) C E- {s) and E+ {t) C £+(5). 

Moreover, for t > 0, v{t) — ±1 quasi-everywhere on E^{t) and ^^i^^-^ ^ E^{t), up to a set 
\dYv u(t)\-negligihle. In particular 

dhiu{t)'^{E-{t)Y > 0, divuit)^{E+{t)Y < 0, div u{t) ^ {E+ (t) U E- {t)y = 0. 

Here, for a Radon measure p and a Borel set E, p\_E denotes the measure defined by 
p^E{B) :^ p{EnB). 

Proof. The first two assertions, as already observed, follow from (|17p and the definition of Ei^ . 
We know that the solution of equation (ITT]) satisfies d^u = —d°'D{u{t)) = Vu(i) for any t > 0, 
but the right-derivative of u = uq -f Vw(t) is nothing else as lim;i_j.o V[u;(i + h) — w{t)]/h. We 
easily deduce that v{t) = liniii^Q[w{t + h) — w(t)]/h (which converges in i?g-strong). Since when 
X G E^{t), w{t,x) = t and 'w{t + h,x) = t + h for h small enough, we deduce that v{x) — 1 on 
that set, in the same way u = 1 on E^{t). 

Observe that the Euler-Lagrange equation for (IT^ is the variational inequality 

{uf)+Vw{t)) ■ {tVv -Vw{t))dx > 0, 



for any v £ Ho{A; [—1, 1]). In other words since u{t) = uq + Ww{t), 

u{t)-V^ > - / uit)-Vv 

t J A 

for any |t;| < 1, and we recover that —\i'w{t)/t £ &D{u{t)). 
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Hence (using Lemma IXT|) . 8^^^-^ C E^{t). Now, if u € H^{A; [-1, 1]) with v ^ ±1 on , one 
deduces that for any s > t, 



'Vv-u{s)dx = 'D{u{s)) . 

J A 

Sending s — > i, it follows 

Vv-u{t)dx > V(u{t)), 



L 



A 

i 

(t) 



hence v € d'D{u{t)). We deduce that ^ Ef{t), invoking Lemma [3?T] □ 



Remark 3.4. We might find situations where \v{t)\ = 1 outside of the contact set. For instance, 
assume the problem is radial, divuo is positive in a crown and negative in the center. Then one 
may have that is a crown {w should be less than t at the center) and E^ is empty. In that 
case, V should be equal to one also in the domain surrounded by the set 

We show now another simple comparison lemma: 
Lemma 3.5. Let uq and Uq in i2(A;M^) such that 

div Uq < div Uq 

in H^^{A). Then for any t > 0, w'(t) < w{t), where w'{t) and 1^(1) are the solutions of the 
contact problem (1121) . the first with Uq replaced with Uq. 

Proof. Let t > 0, e > 0, and w"^ be the minimizer of 



. 1 
mm — 

|«i|<t 2 



\7w\^ dx — / ii;(divitQ — e) 
A J A 



which of course is unique. We now show that < 'w{t) a.e., and since w"^ — >■ w'{t) as e — the 
thesis will follow. 

We have by minimality 

^ dx - I w(i)(divuo) < / ' dx - / {wit) \J w^){diY Uff) , 



A 2 



^ dx - w''(divu^-e) < / ' dx - / A u;'')(div - e) , 



A J A JA ^ J A 



where we denote w(t) V :— max{w{t),w'^} and w(t) A := mhi{w{t),w'^}. Summing both 
inequalities we obtain 



{w{t) y w" - w{t)) div Uq < {w" ~ w{t) Aw" ) {div Uq- e) , 
A J A 

from which it follows e /^(w^ — w{t))^ dx < 0, which is our claim. □ 
Corollary 3.6. Under the assumptions of Lemma \3.5[ 

E-{t) C E'^{t) and E'^{t) C E+{t), (21) 
and it follows that v'{t) < v{t), for each t > 0. 
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Proof. Eqn (PT|) follows at once from the inequality w'{t) < w{t) (Lemma 13. 5p . We deduce, of 
course, that also {t) C E'^(t), and E'^{t) C E^(t). Consider the function v — v'{t) Av{t) — 
mm{v'{t),v{t)}. As it is ±1 on E'^{t), it follows from Lemmas 1531 and [XT] that -Vu G dV{u'{t)). 
In the same way, v' — v'{t) V v{t) — max{u'(t), is such that — Vw' G d'D{u{t)). Since 

/ |Vw|^dx+ / iVv'l^da; = / \\/v{t)\^dx+ [ \S/v'{t)fdx, 

J A J A J A J A 

either |Vwp dx < \Vv{t)\'^ dx or |Vw'p da; < |Vw'(i)p da;. By minimality (as -Vv{t) = 
d°V{u{t))) it follows that v ^ v{t) and v' = v'(t). □ 

3.2 The support of the measure diwu 

Throughout this section we will assume that div uq is a bounded Radon measure on A. 

Lemma 3.7. Let uq e L'^{A;R'^) n domi?, d > and u ^ {I + (552?) "^(uo). Then for a positive 
Radon measure ^ G H^^{A), the Radon- Nikodym derivatives of div u and divwQ with respect to 
fi satisfy (div u//i)(a;) < (divMo/^)(a;) for /i-a.e. x G and (divu//i)(a;) > {divuo/ ^){x) for 
fi-a.e. X G £^ . In particular, divu << divuo and (div it)* < (divuo)*. 

Remark 3.8. It follows from the Lemma that divu = 6div uqI— U for some weight 
6{x) G [0, 1]. We can build explicit examples where < 1 at some point. Consider for instance, in 
ID, A = (0, 1) and the function ua{x) = if a; < 1/3 and a; > 2/3, and 2 - 3a; if 1/3 < x < 2/3. 
Then, one shows that u{t) is given by 



u(t, x) 



3t if X < i 



1 - 2V3i if i < a; < a{t) ■= ^ + 

2 - 3a; if a{t) <x< b(t) := 1 - 



yT+6i-l iix>b{t) 



until t = 1 — 2^/2/3. We have div u{t) ~ u{t)x = (1 — — 3t)(5i/3 — 3x(a(t).6(t)) for such t: 

^u(t) ~ {^/'^} stays constant for a while (and disappears suddenly right after t = 1 — 2a/2/3), 
while the density of the measure div u{t) goes down monotonically until it reaches zero (notice 
that v{t) will jump right after 1 — 2\/2/3), while E~f^^^ = [ait), P(t)) shrinks in a continuous way, 
and carries the constant continuous part of the initial divergence (—3). 

Proof. We have u — uq + 6\7v with — Vw G d'D{u). Let a; G f^J". Recall that the precise represen- 
tative of V is defined by 

v{x} Imi ^ , 

p-*-0 OJNP 

where oj^ = |i?(0, 1)|, and that this limit exists quasi-everywhere in A. We assume also that 
v{x) = 1. 

Then, for a.e. p > 0, one may write 

divu — u ■ V dTi.^ 

B(x,p) JdB(x,p) 



[ UQ -vdl-L^ +5 I Vv ■ udn^ 

JdB{x,p) JdB{x,p) 

I divuo + 5 ( Vu -vdl-L^. (22) 

Jb{x,p) JdB{x,p) 
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Now, let f{p) — ^) Jg^i^ ) vdH^ (which is well-defined for any p). Then, since v{x) — 1 



ldB{x,p) 

and V < 1 a.e., 



limsup/(p) = Nujn 



One can also show that for a.e. p > 0, f'{p) = [1/ p^^^) /g^^^ ^-j Vu • vdH^, in fact / is locally 
in some small interval (0, po). 

Since v < 1 a.e., /(p) < Nujn a.e., so that 



liminf / , / Vv-vdT-L^dr = liminf / f'(r)dr 

= liminf/(p)-/(e) < 



e-i-O 



for any p. If follows that for any p small, the set = {r G [0,p] : jgsix r)'^''^ ' ^ — 0} 
has positive Lebesgue measure, and for any r e /+, we deduce from ^22\i that Jg^^ ^^divu < 

Now consider p a positive Radon measure: /x-a.e., we know that the limits 

divu, , /B(...)div?/ divwo, , /s(.,r)divMo 

(x) = hm — — -— and (a;) — hm — - — -— 

/X 1-^0 p{B{x,r)) p r^o p{B{x,r)) 

exist. If moreover, as before, x e £^ and v{x) = 1 (which holds p-a.e., since p G H~^{A)), 
we can find a subsequence r„ such that /gj-^, ^ ■) divw < /g^^ ^ divwo for each n, and it follows 
(div u/ p){x) < {div uq/ p){x). □ 



The following corollaries follows: 

(t) - 



Corollary 3.9. Let t > s > 0: then (divu(i))^ < (div 7i(s))^. In particular, ^^^^^ ^ £^ 
|divw(s)|-a.e. in A. 

Proof. Indeed: if i > s, then u{t) = {I + {t — s)d'D)~^ {u{s)). We deduce that for quasi-every x £ 

1 = 6ldiv„(t)(a;) < (divM(s)/(divu(t))+)(a;), and it follows (divM(t))+ < (divu(s)/(divu(t))+)(divu(i))+ < 
(divu(s))+. □ 

Corollary 3.10. We have that (divu(t))* ^ (divito)^ as t — >■ 0, weakly-* in the sense of 
measures. Moreover, £^^^ C £'^(0) (up to a \div uo\-negligible set), and divuo I— (-E'^(O)) > 0, 
divMoL(i;-(0)) < 0. 

Proof. We know that as i — ?> 0, u(t) — >■ uq in L'^{A; K^), and thanks to the boundedness of div u{t) 
it follows that div u(t) divito in the sense of measures. Now consider a subsequence (tk) such 
that (divu(tfe))+ p, (div u{tk))^ — ^ v. Since p — v = divuo, it follows that p > (divito)^ and 
v > (div Mo) ^. The reverse inequalities follow from Lemma 13.71 and the first part of the thesis 
follows. 

From the previous results we obtain that for each t, one can write 

(divM(t))+ = e't(x)(divuo)+ 

The function 6t{x) ~ liminfp_>.o(/g^^ divM(t)+)/(Jg^^ ^^(divup)^) is well-defined on the set 
which supports the measure (divuo)^, and we find that Ot{x) < 1 is nonincreasing in t. 
Hence there exists for all x e the limit limf_>.o 0t{x) ~ supj^Q 9t{x), and this limit must be 1 
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(div uo)^-a-e-, otherwise this would contradict that (divw(t))"*' ^ (divwo)'*'. It foUows that up to 
a (divMo)+-neghgible set, C Ut>o{2; e ' ^*(^) > 0}- 
Now, if a; e and 9t{x) > 0, then x £ ^t(ty indeed, 

/b(.,p) div^W ^ (div^ft))+(g(x,p)) - (divuft))-(g(x,p)) p^o 

/^(^^^)|div«(t)| (divZi(i))+(i?(.T,p)) + (div«(t))-(i?(x,p)) ' 

since 

{d\Yu{t))-{B{x,p)) < {diY u-){B{x,p)) 

= o((divw+)(i?(x,p))) < o{{divu{t))+{B{x,p))) 
(the equahty is because x G S^^, the last inequality because Ot{x) > 0). It follows that 



and the conclusion follows from Lemma [3.31 □ 
3.3 The regular case 

Let us now assume that divuo — g £ LP{A), p > 1. The obstacle problem which is solved by w{t) 
can be written ^ 

min — / \\/w{x)\'^ dx — / g{x)w{x)dx. 

weH},:\w\<t2 Jji J A 

Standard results show that w{t) e W^^^p(A), (see Theorem 9.9 in [10]). In particular, we have 
that in the sense, 

-Aw(t) = 9X{\w(t)\<t} 
and, since u{t) = uq + \'w{t), we deduce that in this case 

div u{t) = divuoXE{t) (23) 
for any t > 0. In particular, formally, we deduce from (|lip that 

div wo = Av{t), (24) 

and since Av{t) is the jump of the normal derivative of v{t) on dE^{t), we find that these sets 
shrink with a normal speed |Vw(i)|/|divuo|. 

This can be written rigorously in the sense of distributions: {E^,E^,v) are such that v G 
L^{[0,T);H^{A; [-1,1])), v = ±1 on for a.e. t and x, and for any (j) e C^{[0,T) x A), 



drvuo{x)(l){0,x) dx + / / divito(a;)x£;(t) (a;)-^(a;, i) dx 
A Jo J A 

/ Vv{t,x) ■V(l}{t,x)dxdt = 0. (25) 

/O J A 

We observe that the evolution equation (P5|) is reminiscent of the enthalpy formulation of the 
one-phase Stefan problem ]TE\ . 

We expect that with either the additional information that divuo is a.e. nonnegative on E'^ 
and nonpositive on E~ , or that the maps E^{t) are nonincreasing, then (|25p characterizes the 
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unique evolution (jlip . On the other hand, without this additional assumption, then a time- 
reversed evolution with will satisfy the same weak equation, with uq replaced with —uq. With 
both assumptions we can actually show the following result: 

Proposition 3.11. Let E+,E^ be measurable subsets of A x [0,T], and v € L'^{[0,T)] H^{A)) 
with \v\ < 1 a.e., v = ±1 a.e. on E^ , and satisfying (j25p . Assume in addition that idivug > 
a.e. on E^ , and 

E^{t) C E^{s) fora.e.t>s. (26) 
Then u{t, x) :~ Uq{x) + V w(s, x) ds is the unique solution of 

Proof. Let w{t) = J^v{s)ds. Thanks to (pS)) . we have that < t for a.e. x & A, and 

w{t,x) = ztt for a.e. x g E^{t), for all t. We can approach test functions of the form X[o,t]4>{x)y 
(f) E Hq{A), with smooth functions and pass to the limit to check that 

/ divuQ(j)dx— / divuQ(j>dx = / Ww{t) ■ S/cjjdx , 
J A JEit) J A 

for almost all t (up to a negligible set, which we can actually choose independently of </>, as Hl{A) 
is separable). 

If we choose 4> — w{t, •) as the test function in this equation, we find 

div UQ{x)(f>{x) dx ~ / div UQ{x)w{x,t) dx — / div Uo{x){(f>{x) — w{x,t j) dx 
A J A J E{t) 

Ww{t,x) ■W4>{x)dx - / \Ww{t,x)\'^ dx , 

A J A 

^ |Vw(t,a;) - V0(a;)pdx+ i / \\7 (f>{x)\'^ dx - ^ [ \\7w{t, x)\'^ dx . 



Ja Ja J a 

If |0| < t, we have that — divuo(x)(0(x) — w{x,t)) > for a.e. x E E{t), so that w{t) is the 
minimizcr of (fT2l) and the thesis follows. □ 



Remark 3.12. As mentioned above, it is a natural question whether assumption (126)1 is necessary 
to prove this result. For instance, in case E^ and E~ are closed sets in [0,T) x A with E^{t) n 
E~{t) = for any t > 0, and {divuo = 0} is a negligible set, then one can actually deduce (|26|) 
from (j25p . Indeed, using localized test functions 0(a:)x[s,t]) one shows first that v is harmonic in 
A \ E{t) for a.e. t, and then that Se(s) divuo'/'c^a^ — /g^j-j d\vuo4>dx > 0, and ([26]) follows. 

Remark 3.13. When p > N/2, we can deduce some further properties of w from the regularity 
theory for the obstacle problem [7]. Indeed, letting * E Hg{A) n W'^'P{A) such that -A* = g, 
we have that w — w — E Hq {A) solves the obstacle problem 



min - f \\/w(x)\'^ dx. 



-t-*<i5<t-* 2 

Since p > N/2, we have w{t) E C"(A), with a = 2 - N/p, so that E{t) = {\w{t)\ ^ t} is a 
closed set. In this case, v{t) can be defined as the harmonic function in A \ E(t) with Dirichlet 
boundary condition v(t) — on dA and v{t) — ±1 on E^{t). Moreover, it is easy to check that 
— Vw(t) E d^'D{u{t)), and v{t) is continuous out of the singular points of dA U dE{t). 
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Remark 3.14. li A = one can easily show easily by a translation argument that uq G 
H^{A;'R'^) u{t) G H^{A;R^) with same norm, so that the 77^-norm of u{t) is nonincreasing. 
In this case, ^^(^j is a. e. -equivalent to the support of (divw)^ and since from the equation it 
follows u = uq a.e. on (since v — ±1 a.e. on , so that Vw = a.e., the problem being 
in general that this will not be true quasi-every where), we deduce that divw = divuo a.e. on 
E+UE^ ^ spt(divu). 

4 Examples 

4.1 The antiplane case in dimension 2 

Let = 2 and fc = 1. We have 

J(V') = |rot^|(A) = supj^V^-V' : «eCr(A;[-l,l])| 

where rot?/; = 9i?/'2 — d2ipi and = {82,— di). Then, we check easily that in L'^{A;M.'^) the 
functional J is the support function of the closed convex set 

K = {V^« : veH'oiA;[-iA])}. 

As we mentioned in the Introduction, this functional appears as limit of the Ginzburg-Landau 
model in a suitable energy regime [19] . 

Letting ■0^ — {ip2, —i'l), we get J{ip) = |div?/;-'-|, so that the flow can be described as above. 

Proposition 4.1. Let uq G L^{A;'R^) with rotuo = .9 G LP (A), p > 1. Then for t > there exist 
nonincreasing left- continuous closed (and disjoint) sets E^{t) C {±g > 0}, such that TOtu{t) = 
rotMo(x_E-(t)u£;+(t))- Moreover, letting E^ = Ut>o{i} x E^{t), there exists a function v{t, x) with 

V = ±1 a.e. on E^ such that {E~^,E~,v) are the unique closed sets and function solution of the 
weak Hele-Shaw flow (P5)) . 

4.2 The one-dimensional Total Variation Flow 

Let now = 1, fc = 0: the previous analysis also provides interesting qualitative information on 
the behavior of the flow of the Total Variation, in dimension 1. 

We consider uq G L'^{{a, b)), a < b, and the flow u{t) of the total variation J{u) := sup{/^ uv' dt : 

V G C^{a, b; [—1, 1])}. Notice that in this situation, the function w which minimizes ((HI), being in 
Hoia, b), is also in C^/^([a, b]) with w{a) — w{b) — 0. In particular, the sets E'^lt) defined in (fTS]) 
are closed, disjoint sets compactly contained in (a, 6). 

We can state the following result. 

Proposition 4.2. The function u{t) is the unique minimizer of 

minJ(u) + — \u — uo\'^dx. 

Moreover there exist nonincreasing, disjoint closed sets E^{t) C (a,b) such that u{t) — Uq a.e. on 
E^{t), Uq is nondecreasing on any interval contained in E^{t), nonincreasing on any interval 
contained in E^{t), and u{t) is constant on each connected component of (a, b) \ {E^{t) U E^{t)). 
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If uq is smooth enough, one can also characterize the speed of the boundary points of (t) 
in term of uq and the size of the intervals of (a, b) \ U E^{t)). 

Proof. The first part of the thesis is a consequence of Remark 13.21 Then, if uq G BV{a,b), the 
thesis is a consequence of Lemma [3.31 Indeed, for a.e. x on E^{t), we have dxw{t,x) = and 
u{t,x) — uq{x) + dxw{t,x) — Uo{x). If / C E^(t) is an interval, since the measure Du{t)\_I 
must be nonnegative, u{t) is nondecreasing on /, but as u{t) = uq a.e. on I it follows that uq is 
nondecreasing on /. 

If Uq ^ BV{a,b), we use the fact that for all e > 0, u{e) E BV{a,b). Then the Proposition 
holds for t > £, and we have u{t) — u{e) a.e. on E^{t), u(e) is nondecreasing on any interval 
contained in E~^(t), nonincreasing on any interval contained in E~{t), and u{t) is constant on 
each connected component of (a, b) \ {E~^{t) U E~{t)). The sets do not depend on e, as they are 
defined as the contact sets in (fT2)) . Sending then e — we deduce the result. □ 

We can deduce the following, quite interesting result — see also [17[ [3 [14] for other results on 
the one-dimensional Total Variation flow and in particular [17, Prop. 4] for a similar statement. 

Corollary 4.3. Let uo = uq + n where uq G BV{a, b) and n is a stochastic process (a, b) with 
n G Lp'{a,b) a.s. and such that \Dn\{I) = +oo for any interval I C (a, 6), almost surely. Let u{t) 
be the total variation flow starting from uq. Then almost surely, att > 0, there is "staircaising" 
everywhere in the interval (a, b): u{t) is constant on each connected component of an open set A{t) 
which is dense in (a, 6). 

Remark 4.4. The property that \Dn\{I) = +oo for any interval /, almost surely, is satisfied 
for instance by the Wiener process (as its quadratic variation is positive a.s.). For a Gaussian 
stationary process, it will depend on the behaviour of the autocorrelation function and can be 
characterized by conditions on the power spectrum of the process, see for instance [4^ for (non 
sharp) conditions. 

Proof. We let A{t) = (a, 6) \ {E^{t) U E^{t)), and from the previous result we know that u{t) is 
constant on each connected component of A{t) while u = uq on (a, b) \ A{t). Now assume there is 
an interval / with / fl A{t) = 0: without loss of generality we may assume that / C E~^{t). Then 
Mo must be nondecreasing on /, in particular there exists L' <Z I with |Z?mo|(/') < +oo. But this 
yields that \Dn\{L') < +00, which is a.s. impossible. □ 
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